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Abstract 

Ab initio wave-function-based methods are employed for the study of quasiparticle energy bands 
of zinc-blende ZnS, with focus on the Zn 3c? "semicore" states. The relative energies of these states 
with respect to the top of the S 3p valence bands appear to be poorly described as compared to 
experimental values not only within the local density approximation (LDA), but also when many- 
body corrections within the GW approximation are applied to the LDA or LDA+U mean-field 
solutions [T. Miyake, P. Zhang, M. L. Cohen, and S. G. Louie, Phys. Rev. B 74, 245213 (2006)]. 
In the present study, we show that for the accurate description of the Zn 3d states a correlation 
treatment based on wave function methods is needed. Our study rests on a local Hamiltonian 
approach which rigorously describes the short-range polarization and charge redistribution effects 
around an extra hole or electron placed into the valence respective conduction bands of semicon- 
ductors and insulators. The method also facilitates the computation of electron correlation effects 
beyond relaxation and polarization. The electron correlation treatment is performed on finite 
clusters cut off the infinite system. The formalism makes use of localized Wannier functions and 
embedding potentials derived explicitly from prior periodic Hartree-Fock calculations. The on-site 
and nearest-neighbor charge relaxation lead to corrections of several eV to the Hartree-Fock band 
energies and gap. Corrections due to long-range polarization are of the order of 1.0 eV. The dis- 
persion of the Hartree-Fock bands is only little affected by electron correlations. We find the Zn 
3d "semicore" states to lie about 9.0 eV below the top of the S 3p valence bands, in very good 
agreement with values from valence-band x-ray photoemission. 
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I. INTRODUCTION 



The accurate ab initio treatment of many-body effects in the ground and excited states 
of sohds remains a challenging problem in physics. For the ground-state (GS) properties 
of periodic systems reliable ab initio methods are at present available. Two different for- 
malisms can be distinguished: approaches based on Kohn-Sham (KS) density functional 
theory (DFT)^ and methods which rest on explicitly computing many-body wave functions 
(WF's)^'^, by means of quantum chemical techniques^. The accuracy achieved by the two 
types of methods for the GS properties is often comparable. For the case of energy bands 
and excited states, however, the results obtained by KS-DFT and WF-based approaches 
are different. Although DFT formally is a GS theory, it has been commonly employed in 
the computation of energy bands, i.e., excited states, and in some cases yields surprisingly 
good results, e.g., when hybrid density functionals or exact exchange schemes are utilized. 
Nevertheless, a rigorous justification for assigning the KS eigenvalues or KS energy differ- 
ences to excitation or quasiparticle energies of correlated electron systems is still to be given. 
Since there is no solid basis for such an extension of GS KS-DFT, when inaccurate or wrong 
excitation energies are obtained, they can not be improved in a systematic and controlled 
manner. 

The widely used local density (LDA)^ and generalized gradient approximations (GGA)- 
to the KS-DFT provide a rather poor description for binding energies and band gaps of 
semiconductors and insulators. The self-interaction error introduced with the LDA or GGA 
potential and the absence of discontinuity in the exchange potential have been identified as 
the main reasons for that. The incomplete cancellation of the self-interaction in LDA or 
GGA is also regarded as a major cause for the underestimation of the binding energies of 
"semi core" d states in systems such as the 11'^- VI"^ semiconductors''"^. An accurate ab initio 
description of the electronic structure of these seemingly conventional sp semiconductors 
has appeared to be a rather difficult task. 

Numerous refinements of KS-DFT have been suggested to enable the treatment of many- 
body effects in correlated electron systems. Many-body corrections to the KS mean-field 
eigenvalues as obtained within the GW approximation^*^ to the self-energy often yield good 
agreement between theory and experiment with respect to the band gaps of various group- IV 
and III-V semiconductors. The GW approximation is established as the standard method for 
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calculating quasiparticle energy bands in those materials. It does fail, however, to describe 
complex many-body effects in d-electron systems such as the II^-VI^ compounds^iii"— . The 
quasiparticle energies of the d states, for example, are about 2.0 eV higher than the exper- 
imental valuesi^ii^. This result remains the same when the mean-field LDA+ [/^^ Hamilto- 
nian is utilized^ as a basis for the GW treatment-'^. The problems were attributed to the 
self-interaction errors, still present in the standard GW approximation^^. 

It has been argued that determining self-consistently the optimum mean-field Hamiltonian 
rather than utilizing the LDA {+U) would largely improve the description of quasiparticles 
within the GW approximation^^. Indeed, the self-consistent GW method of Kotani et 
ai— yields in the 11^- VI"^ and III^-V"^ compounds d-state quasiparticle energies in good 
agreement with experimental values, i.e., deviations of at most 0.8 eV are observed. Yet, 
the approach overestimates systematically, although by little, the band gaps. This was 
attributed to the lack of excitonic effects in the dielectric function. 

Methods which treat the exchange energy within the exact exchange (EXX) scheme and 
are thus considered to be self-interaction free have been developed as well. For example, the 
full-potential EXX KS-DFT calculations of Sharma et al.^ including core- valence exchange 
interactions yield very good agreement with the experiment for the binding energies of d 
electrons in ZnS. The band gaps of the II^-VI"^ compounds and other sp semiconductors, 
such as Si, Ge, and GaAs, are, however, largely overestimated in that computational scheme. 
This trend is reversed in insulating compounds such as BN, C and some nobel gas solids, 
where the gaps are underestimated^. Excluding the core-valence exchange interaction, on 
the other hand, brings the band gaps close to experimental estimates, but it simultaneously 
worsens the binding energies. EXX calculations within the pseudopotential framework yield 
band gaps and cohesive properties of various sp semiconductors in very good agreement with 
the experiment^. The binding energies of the d states, however, remain wrong in such a 
treatment^. 

A different approach to energy bands in periodic systems is based on ab initio WF meth- 
ods. Their major virtue is the use of well-defined and controllable approximations. Con- 
verged results for binding energies, lattice parameters, and band gaps can be here achieved 
by systematic improvement of the many-body wave function^. The disadvantage of those 
formalisms is the larger amount of work involved as compared to the DFT-based calcula- 
tions. The large computational effort is related to the physics of the problem. When an 
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electron is added to a conduction-band or removed from a valence-band state of a semicon- 
ductor or insulator, the correlation hole of the extra particle (hole) is long-ranged because 
of polarization and relaxation effects. The corresponding quasiparticle, i.e., the particle plus 
its correlation hole, is moving in form of a Bloch wave through the system. The computation 
of the energy dispersion associated with the quasiparticle requires an accurate description of 
the correlation hole. The latter differs considerably from the correlation hole of an electron 
in the ground state for which the correlations are of van der Waals type and do not involve 
long-range (LR) polarizations. The different nature of the correlation holes in the ground- 
and (A^±l) states is simply ignored in the DFT methods and hence the large computational 
effort is avoided. Taking into account this difference is essential for the correct determination 
of energy bands and gaps. 

In contrast to ground-state calculations^Si"— , energy band calculations based on quantum 
chemical methods are still limited to non-metalic systems^i^iS^i^Sr— . They rest either on 
second-order M0ller-Plesset theory (MP2)^""— , coupled-cluster (CC) techniques^ or effective 
local Hamiltonian (LH)^"— i^i^ approaches. Some of these formalisms make use of the local 
character of the short-range part of the correlation hole^^'^^'^'^'^""— , utilizing thereby local 
operators^'^ related to sets of real-space, localized orbitals. The localized orbitals are either 
derived by means of various localization procedures^^^^ applied to the canonical Bloch 
orbitals of the self-consistent-field (SCF) GS wave function^^'^'^'^"— or obtained directly 
from a crystal Hartree-Fock— calculation carried out in the Wannier representation^^'^'^. 
The local orbital schemes of Pulay and Saeb0 utilized in local M0ller-Plesset and coupled- 
electron-pair correlation approaches to finite systems^i^ have also been adopted for deriving 
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local orbitals for crystalline compounds, see, e.g., Ref. 

In the present study, we compute the correlation-induced corrections to the HF valence 
and conduction energy bands of c-ZnS, a representative of the the 11"^- VI"^ compounds. There 
is much interest toward the electronic structure and properties of the 11-^- VI"^ systems due to 
their optical properties^"— and potential exploitation in optoelectronic devices and medical 
instrumentation. For example, thin films of these materials grown on substrates of GaAs 
give rise to interfaces emitting intensive white light, i.e., semiconductor-based light emitting 
devices^i^. Moreover, some properties of mixed phases of the type ZnS^^Sei-a; such as the 
width of the fundamental band gap may be continuously tuned depending on the explicit x 
valued. 
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The goal of these WF-based investigations is twofold. On one hand, we extend the 
area of applicability of the WF-based band structure calculations, building on previous 
investigations performed for the simpler compounds MgO-^ and BN— . On the other hand, 
our study attempts to provide new insight into the many-body effects which operate in this 
specific 3d electron system. Of particular interest are the binding energies of the "semicore" 
Zn 3d states since their accurate treatment remains out of reach for standard DFT and GW 
methods, see Refs. Q, 11 



and the discussion above. HF calculations carried out within the 
linear-combination-of-atomic orbital method^ or the full-potential linear-muffin-tin-orbitals 
approach^ overestimate as well, by about 3-5 eV, the binding energies of d states in II^-VI"^ 
and III'^-V"^ semiconductors. 

Our WF-based studies are performed in the framework of the LH approacb^^""— and 
quasiparticle approximation^, where the quasiparticle energy bands are expressed in terms 
of matrix elements of an effective local Hamiltonian. A short overview of the method is 
presented in Section II, since a detailed account has already been published^'^Si'^'^'^. 

II. THEORETICAL FRAMEWORK 

1. Quasiparticle energy bands 

The correlated quasiparticle valence- and conduction-band energies of periodic (closed- 
shell) systems are defined qs^-^'—^^ 

e^,. = i?o^-«-;ii^K-;), 

V ' 

e^.. = {^ll'\HKJ.')-E^, (1) 

^ V ' 

where is the energy of the correlated ground state of the neutral A^-electron crystal and 
^kiZ^ and E^^^ are the correlated energies of the crystal (A^+1) and (A^— 1) electron states, 
respectively. The quasiparticle states I'^^J'J) of the (A^+1) electron system are related to 
the corresponding SCF counterparts I^^i/ct^) througb^'^'^ 

KV) = ^K^;), (2) 

where u is a band index, and a and k are the spin and momentum of the (A^+1) state. The 
wave operator Q in a. coupled cluster formulation is given hj Q = , where is a scattering 
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operator. A similar relation holds for the (A^ — 1) states of momentum — k, band 

index fi, and spin —a obtained when an electron is removed from the system. 

The SCF (A+1) states are expressed in terms of local electron-addition single-particle 
configurations 



KJ-J) = 4f E «-(k)e^'-^^ <„.l$sc.), (3) 



n,Rr 



where No is the total number of unit cells and a^„(k) is the structure matrix^. The local 
operator wl^^,^^ creates an electron in a real-space localized conduction-band Wannier orbital 
(WO) IwnaCR-i)) of index n, centered at a site defined by the lattice vector R/ in unit cell 
/. The SCF GS wave function |$scf) and energy E^'^'^ of the A^-electron system as well as 
the WO's are derived via a periodic HP calculation for the extended system. Analogously 
to the HF (A^ + 1) states, the HP states are expressed in terms of local single-hole 

configurations I^R^mo-) ^y utilizing local annihilation operators WR^mo- for the valence-band 
WO's \wma(R-i)) ■ By combining Eqs. ([2]) and ([3]) with Eq. ([1]), the correlated conduction- 
band energies are expressed in terms of the matrix elements of the local effective Hamiltonian 
H-ff 

R/ nn' 

SCF/ 



X 



H 



— Eq SjijoSnn', (4) 



On'cr,R.jricr 



where Eq = (^scfI^o-^^^oI'^scf) and the wave operator Qq transforms $scf into the true 
ground state \1'^ of the A-particle system. In Eq. HI the matrix elements of H'^^-^ are 
defined in terms of the real-space correlated (A + 1) electron states |^R^V)=^'"^R7ncrl'^scF)- 
A Wannier-like transformation relates I'^^'^na) "with their k-space counterparts I^I/^j,^^)- 

To determine the correlation corrections to the SCF band energies, ^^ua—^^vo-^ 
e^J^, the SCF energies e^J^ are first computed by utilizing the local matrix elements 
HonlRin={^scF\won'cTHwl^^^J^scF)-E^'-''5B,j05„^,- The expressions for the valence-band en- 
ergies, corresponding matrix elements and correlated wave functions are analogous. 
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Because of the predominantly local character of the correlation hole, it suffices to compute 
all relevant matrix elements -f^on^R/n ^om' Rjm using sufficiently large clusters cut off the 
infinite systent^^'^'^i^'^. These matrix elements can be determined by utilizing CC or mul- 
ticonfiguration (MC) and multireference (MR) configuration-interaction (CI) calculations^ 
when correlations are strong. In principle, the correlation treatment can be facilitated by the 
method of increments, where the total correlation energy is assembled from contributions 
originating from selected groups of localized orbitals^iiSli^. 

The scattering operator S is given by a selected set of excitation operators, see, e. g., 
Refs. For relatively weakly correlated systems, a reasonable ansatz for S includes 

one- and two-particle excitations which generate the correlation hole of the added parti- 
cle, i.e., the relaxation and polarization cloud, and account for the loss of ground-state 
correlation8^i^i^i2£. 



2. Computation of the matrix elements of IF^f 



SCF 

R/n' 



To obtain the correlation corrections to the matrix elements of Yl^f^, -^^on'.R/n ~ H^^, 
we calculate the short- and long-range part of the correlation hole of an added electron or 
hole. The short-range (SR) part of the correlation hole consists of intra-atomic and SR inter- 
atomic relaxation and polarization. This response of the nearest surroundings is determined 
by means of SCF orbital optimizations for a finite cluster around the extra electron or hole 
set into a conduction-band or valence-band WO. The WO hosting the added particle is 
kept frozen in the SCF calculation^ for the (A^=f1) excited state, which is referred to as 
the frozen local hole (electron) approximation^^*^. In first-order perturbation theory, the 
SCF orbital relaxation corresponds to single-particle excitations around the frozen WO, 
see, Ref. |67|. The corresponding wave functions I^R^mo-) ^'^d I^r^V)^ t\ms given by 
the optimized SCF states I^R^mo-) I'^r^V)) respectively. For the case of weakly and 
moderately correlated systems, the SCF calculation for the finite region suffices to describe 
the SR correlation hole. Because the local operators are associated with a set of spin orbitals, 
while, in practice, a spatial orbital set is utilized, the spin degree of freedom is accounted 
for by further optimizing the WO with the extra particle within the previously relaxed 



environment, see Ref. 
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The wave functions expressed in terms of relaxed orbitals within 



the nearby surroundings, including the reoptimized singly-occupied WO, are denoted as 
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l^^L) and ll-KJ- 

In distinction to the correlation hole of an electron in the GS of a system, which extends 
only over few lattice sites^'^'^, the correlation hole of an extra particle in the (A^±l) system 
includes a substantial long-range (LR) polarization cloud. This effect originates from the 
polarization of the farther surroundings and can not be described by calculations on finite 
clusters. Yet, only the SR part of the correlation hole is essential when the energies of 
different (A^+1) or (A^— 1) states are compared because the contribution of the LR polarization 
cloud cancels out. The LR polarization cloud is treated here within a dielectric continuum 
approximation using the experimental dielectric constant^i^^i^^i^i^. Its contribution to the 
correlation-induced corrections to the matrix elements of if^-^-^ is explicitly incorporated and 
hence ionization potentials and electron affinities can be evaluated as well. 

Correlation effects beyond charge relaxation and polarization, including the loss of 
ground-state correlations, are computed by means of single and double CI (SDCI) 
calculations^ on top of the SCF orbital relaxation, separately for each l^^J^a) and 
or l^^'^ma) and state. The loss of ground-state correlations is a measure of correla- 

tion contributions that are present in the A^-electron ground state but absent in the (A^ ± 1) 
states, since in the latter some excitations are blocked. To account for such effects, the 
conduction-band WO's to which an extra electron is attached and the valence-band WO's 
out of which an electron is removed are always kept singly occupied in the SDCI wave func- 
tions. To treat in a unified manner both valence- and conduction-band states and to avoid 
spurious charge flows within the finite cluster C, the infinite frozen crystal environment of 
C is explicitly incorporated in the calculations. The effect of the electrons in the localized 
occupied orbitals of the environment £ on the electrons in the loca. 
sumed in an effective one-electron potential V'^"^'', see, e.g., Refs. 
potential V^"^^ is computed by utilizing the crystal self-consistent Fock operator F'^^v^^ and 
the Fock operator of C, F[P'^], associated with the density operator pcj5^^_ ^^le cluster 
calculations, yf^'^f'=F'^ry'^t _ F^p*^] is added back to F[P''] and the correlation calculations 
are thus effectively performed in an infinite frozen HF environment. 

Our embedded cluster formalism requires Wannier-Boys localization^® of the HF core, 
valence, and lowest-lying conduction-band states. The optimally localized valence-band and 
conduction-band WO's which are centered at lattice sites within C are then projected onto 
the set of Gaussian-type-orbital (GTO) basis functions assigned to those cluster sites^^i^i^. 



ized orbitals of C is sub- 



70, lDJ- The embedding 
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Their longer-range tails that involve Gaussian basis functions at sites outside C and or- 
thogonal to the basis functions in C are thus cut off. Hence, the projected WO's that are 
centered at sites within the central part of C bear closest resemblance to the original WO's, 
while those near the edges of C are represented asymmetrically in the basis set of C. For 
the post-HF treatment, the cluster C is thus subdivided into a central active region Ca, and 
a buffer region Cb- The role of the sites and corresponding basis functions within Cb is 
to ensure an accurate representation of the longer-range tails of the WO's centered within 
Ca- For sufficiently large Cb regions, the original WO's in Ca remain in practice unaffected 
by the projection, i.e., their projected counterparts show the same properties as the origi- 
nal ones^i^*^. In the post-HF calculation, the occupied and lowest-lying conduction-band 
WO's centered within Ca are explicitly correlated while the WO's within Cb are held frozen. 

The projected WO's are neither normalized nor orthogonal. Thus, they are groupwise 
orthonormalized in the order active core, active valence, active lowest-lying conduction- 
band orbital groups, buffer core, buffer valence, and buffer lowest-lying conduction-band 
orbit als^i^i^. Our localized (A^± 1) states are expressed in terms of those projected WO's. 

To generate the virtual orbital space, the projected atomic orbital (PAO) formalism of 
Pulay^^ and further elaborated by Hampel and Werner''^ is adopted. The virtual PAO's are 
constructed from the GTO basis functions assigned to the sites in Ca after projecting out 
the orthonormalized occupied and lowe st-ly ing conduction-band WO's of C via a Schmidt 
orthogonalization, see also Refs. |6l| and l62l They are subsequently Lowdin orthonormalized 
among themselves. This procedure yields thus PAO's orthogonal to those occupied and 
lowest-lying conduction-band WO's. While the PAO's in Ca are identical with the virtual 
orbitals of the crystal, the PAO's in Cb deviate from those of the infinite system. In the 
correlation calculations, however, the variational virtual orbital space consists only of PAO's 
within Ca- 



III. CRYSTAL STRUCTURE OF ZINC-BLENDE ZNS AND COMPUTATIONAL 
INFORMATION 

ZnS is representative of the 11-^- VI"^ semiconductors crystallizing at ambient pressure in 
the cubic zinc-blende structure. The zinc-blende form is face-centered (fee), with the F43m 
space group symmetry. The corresponding lattice constant is 5.415 A''^. The local point- 
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group symmetry for both Zn and S is T(i, with a four site nearest-neighbor coordination. 

GTO basis sets of triple-zeta quahty augmented with polarization functions are apphed 
in our study for the S atoms. We employed the 86-3 IIG* GTO basis set of Lichanot et alJ- 
derived for crystalline MgS but reoptimized in periodic HF calculations the exponents of the 
outermost two sp shells as well as that of the single-Gaussian d polarization function. The 
exponents of the reoptimized sp shells are 0.27202 and 0.125. The reoptimized d polarization 
function has an exponent of 0.435. The latter value is very close to the one of 0.479 used in 
correlated molecular calculations^^. 

The basis set used for Zn is a 86-41 lld3 11 GTO basis set which is derived from the 86- 
4111d41G basis set constructed for ZnO by Jaffe et alM and recently employed for ZnSa;Sei_3; 
by Homann et air- . The basis set of Jaffe et al. was properly modified to achieve a triple- 
zeta basis for the description of the d states. For this purpose, the most diffuse primitive of 
the 3(i shell was decontracted and the exponents and contraction coefficients of the valence 
3(i shell as well as the exponents of the outermost single- Gaussian polarization Asp, 5sp, 
6sp and 4d, bd functions were reoptimized. Their optimized values are presented in Table 
IVl The exponents of the outer sp and d shells underwent almost no variations during the 
HF optimization as compared to those reported by Jaffe et al.. For the Zn 3d states, this 
outcome is due to their "semicore" nature. 

The periodic HF calculations were carried out with the crystal program^. The 
Wannier-Boys orbital localization module in crystal was employed to construct localized 
Wannier functions^*. The projected WO's for the finite cluster C and the matrix repre- 
sentation F^p^^ of the crystal Fock operator in the basis set of C are both obtained with 
the CRYSTAL-MOLPRO interface program^^. In addition, the CRYSTAL-MOLPRO interface 
generates the virtual PAO's. The Fock operator of C, F[P^], and the density operator P'' 
are calculated with the molpro^ suite of programs. The embedding potential for C is also 
computed with the help of molpro. The post-HF correlation calculations are carried out 
with the MOLPRO package as well. 
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IV. CORRELATION EFFECTS ON THE BAND STRUCTURE OF ZINC- 
BLENDE ZNS 



A. HF band structure 

The HF band structure of ZnS, computed with the triple-zeta GTO basis sets described 
in Sec. Ill, is shown in Fig. [1] The upper valence bands have dominant S 3p character with 




r X w K r L w 
k 

FIG. 1: Hartree-Fock band structure of ZnS. 



very small contributions from the Zn 4s and 4p states. The projected density of states (DOS) 
indicate that the lowest conduction bands have predominantly Zn 4s and 4p character, see 
Fig. |2l The four lowest-lying Zn 4s and Ap conduction bands are entangled with a somewhat 
higher-lying dispersive band which has predominantly S 4s, 4p character, see Fig. [1] The 
admixture of S 4s and ip states is apparent in the DOS plots in Fig. O 

The fundamental band gap of the system is direct, T^^ — )■ F^, in agreement with 
experiment^. As expected, the HF result, 10.96 eV, strongly overestimates the actual 
experimental value of 3.83eV— . The average energy position of the d band complex at 
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FIG. 2: Total and projected density of states (DOS) of ZnS around the Fermi level. 

the r point relative to the top of the valence band is -13.3 eV. At the HF level this quantity 
is overestimated by about 4.0 eV as compared to experimental values of -8.97 eV and -9.03 
eV deducted from valence-band x-ray photoemission spectroscopy—. 

As explained above, the orbital basis set employed for our post-HF correlation treatment 
consists of a set of optimally localized WO's. The Wannier-Boys localization^ is carried out 
separately for the core and valence, and lowest-lying conduction bands. Because of the finite 
energy gap between the lower-lying conduction bands and the higher virtual Bloch states, 
the Wannier-Boys transformation for the former states is straightforward. 

The Wannier-Boys transformation for the four lower- lying Zn 4s, Ap conduction bands 
yields a set of four sp^-like hybrids each oriented along one of the Zn-S segments of a ZnS4 
tetrahedron. Each of those hybrids, centered at the Zn site of the ZnS4 unit, has significant 
weight at the nearest-neighbour (NN) S sites and strong antibonding character with the S 
3s and S 3p orbitals. The larger lobe of each of the sp^-like hybrids is tetrahedrally shaped 
in the vicinity of the Zn site. The WO associated with the higher-lying conduction band 
of predominantly S 4s, 4p character is delocalized over four S sites. Our study is restricted 
here to the correlation-induced corrections for the lower-lying Zn 4s and Ap states. Electron 
correlation effects for the S 4s and 4p bands are not addressed in the following. 
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In order to obtain a set of Cartesian type s and p functions, a Pipek-Mezey localization 
procedure^ is applied to the projected Zn sp^-like hybrids. The Pipek-Mezey transformation 




FIG. 3: Zn 4p-like conduction-band WO's after projection onto a [Zni3S28] cluster. 

yields a set of three degenerate orbitals, each oriented along one of the C2 axes of the ZnS4 
tetrahedron and a fourth, lower-energy WO. These projected WO's, derived from the Pipek- 
Mezey transformation, are plotted in Figs. |3] and HI The degenerate WO's resemble three 
Cartesian p functions while the fourth WO is strongly deformed as compared to a s function, 
i.e., it is tetrahedrally shaped in the vicinity of the Zn site. All these four Wannier functions 
have substantial weight at the NN S sites. Whereas the Zn 4p-like and Zn 4s-like WO's 
are rather extended objects, the WO's associated with the sulphur 3s and 3p valence bands 
are well localized around the S sites and the small tails at the NN Zn atoms are hardly 
visible in the plot shown in Fig. [51 The projected Zn 3s and 3p WO's obtained through 
the Pipek-Mezey transformation are slightly deformed Cartesian-like 3s, Sp^, 3py, and 3^^ 
functions. 

Finally, the projected WO's associated with the Zn 3d bands are very compact and bear 
complete resemblance to the t2 {dxy, dxz, dyz) and e {dx2^y2, dz^) orbitals of a tetrahedrally 
coordinated Zn ion. 

The norms of the projected valence-band S 3s, 3p and of the semicore Zn 3d WO's, 
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FIG. 4: Zn 4s-like conduction-band WO after projection onto a [Zni3S28] cluster. 



centered within the Ca regions of the embedded clusters, are always larger than 0.99 of the 
initial HF WO's. The norms of the projected conduction-band WO's associated with the 
Zn 4s and 4p states are always larger than 0.95. 

B. Correlation-induced corrections to the diagonal matrix elements of the local 
Hamiltonian 

Various incremental contributions to the short-range correlation-induced corrections to 
the diagonal matrix elements of the local Hamiltonian -/?o<R/n obtained by using dif- 
ferent, properly designed clusters. Each embedded cluster C is designed such that a ground- 
state HF calculation which includes in the SCF optimization only the projected WO's in 
the Ca region leads to changes in the total HF ground-state energy of less than few tenths of 
meV. This demonstrates that the projected WO's centered within the Ca region are practi- 
cally identical to the original crystal WO's and that the embedding potentials were properly 
constructed. In other words, the HF description of those finite clusters is indeed equivalent 
to the HF description of the periodic crystal. 
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FIG. 5: S 3p and S 3s valence-band WO's after projection onto a [Si3Zn28] cluster. 



1. Short-range relaxation and polarization effects on the "semicore" Zn 3d band states 

In order to compute the effect of short-range correlations on the diagonal matrix elements 
associated with the Zn 3d states, we designed a [Zni3S28] cluster. We are here concerned 
with on-site relaxation plus relaxation and polarization related to the first (four NN S sites) 
and second (twelve next-nearest-neighbor (NNN) Zn ions) coordination shells around the 
Zn site which accommodates a 3d hole. The Ca region, [Zni3S4], of the [Zni3S28] cluster 
contains a central Zn site, depicted as a large black sphere in Fig. [6l the four NN S atoms, 
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pictured as large white spheres, and the shell of 12 NNN Zn atoms, drawn as average-size 
black spheres. The Cb region includes all S nearest neighbors of the 12 outer Zn sites 
within the Ca region. The on-site relaxation and polarization effects for the Zn t2 and e 



Q 




FIG. 6: Illustration of the [Zni3S28] cluster employed for the study of the valence-band Zn 3d hole 
states and conduction-band Zn 4s, 4p electron states. The active region of the cluster consists of 
the [Zni3S4] fragment. The central Zn site in the active region is depicted as the largest black 
sphere. The other twelve Zn atoms in the active zone are shown as average-size black spheres. The 
four S atoms in the active region are drawn as large white spheres, whereas the buffer S atoms are 
shown as smaller white spheres. 

(A^ — 1) electron states amount to -4.04 eV, see the fourth and fifth columns in Table [H 
Those corrections were computed by restricted open-shell HF (ROHF) calculations where 
the doubly occupied orbitals at the central Zn site hosting the hole are allowed to relax 
and polarize in response to the missing charge. The Zn 3d orbital where the hole resides, 
either of t2 or e symmetry, is kept frozen. A modest contribution of about -0.24 eV to the 
overall on-site correlation effects arises from the relaxation and polarization of the valence 
Zn 3s, 3p WO's. As anticipated, the largest effect comes from the relaxation of the doubly 
occupied Zn 3d WO's, that is -3.79 eV. The 3d hole states in first row transition metal (TM) 
compounds are characterized by large atomic relaxation effects due to the core-like spatial 
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TABLE I: Correlation-induced corrections (in eV) to the diagonal Hamiltonian matrix elements 
for the valence-band S 3p, 3p and Zn 3d hole states. Negative corrections indicate an upward shift 
of the valence bands. 







(0) 








S 3s 


S 3p 


Zn 3de 


Zn 3dt2 


On-site orbital relaxation 


-L303 


-0.852 


-4.031 


-4.036 


NN orbital relaxation 


-0.125 


-0.116 


-2.908 


-2.908 


NNN orbital relaxation 


-1.261 


-1.181 


-0.086 


-0.086 


Hole orbital relaxation 


-0.096 


-0.155 


-0.249 


-0.251 


Loss of GS correlation 


1.760 


0.822 


1.455 


1.599 


Long-range polarization 


-1.537 


-1.537 


-1.537 


-1.537 


Total correlation correction^ 


-2.562 


-3.019 


-7.356 


-7.219 



differential correlation effects other than loss of ground-state correlations are not included 

distribution of the 3d orbitals. 

Orbital relaxation and polarization at the NN S sites are also substantial, -2.91 eV, see 
Table HI Relaxation and polarization within the first zinc coordination shell of the central 
Zn ion result into a small additional correction of about -0.09 eV. Finally, the relaxation of 
the singly occupied Zn t2 or e WO where the hole itself resides brings an extra correction of 
-0.25 eV. 

As illustrated in Table HI the short-range relaxation effects are nearly the same for the Zn 
^2 and e {N — V) states. This is due to the very compact nature of the Zn 3d orbitals. This 
compact character of the Zn 3d orbitals is also the reason for a crystal field splitting that 
amounts to only 0.01 eV. Correlation effects beyond relaxation and polarization actually 
enhance the crystal field splitting to about 0.1 eV, see Section IV. B. 3. The small value of 
the Zn t2-e splitting is consistent with the experimental findings, i.e., no splitting of the peak 
corresponding to the Zn 3d semicore states is observed in the x-ray photoelectron spectrum 
of ZnSi^^. 

In the calculations above, we constructed the relaxed ROHF wave functions l^^J^a) 
first approximations to the correlated wave functions starting from the Koopmans 

electron-removal states mo-)- The infiuence of scalar relativistic effects and spin-orbit 
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coupling on the Zn 3d shell is not explicitly considered. The spin-orbit splitting within the 
3d shell of the Zn'^"^ ion is 0.34 eV, much larger than the crystal-field splittings. The effect 
of those spin-orbit couplings, however, is much smaller than the relaxation and polarization 
effects which determine the actual binding energies of the Zn 3d semicore states. Relaxation 
and polarization effects beyond the first zinc coordination shell are evaluated via a dielectric 
continuum approximation, see below. 

The small radial extent of the Zn 3d orbitals and the large separation between the Zn 3d 
and S 3p energy levels leads actually to small mixing between the Zn 3d and S 3p functions. 
Ligand p to metal d charge transfer (CT) effects are therefore expected to be very small 
in ZnS. Trial multiconfiguration, complete active space (CAS) SCF calculations^ for the 3d 
(A^— 1) electron states were also performed, where the central Zn 3d orbitals and the twelve S 
3p orbitals centered within the Ca region were all included in the active orbital space. These 
calculations show indeed that the weight of S 3p to Zn 3d CT configurations in the CASSCF 
wave functions is not larger than 1 %. The energy separation between the highest-lying S 
3p (A — 1) state and the lowest-lying Zn 3d {N — Tj state is also found to be substantial, 
4.5 eV. The optimized active orbitals turn out to constitute a set of optimally localized 
atomic-like active orbitals with small mixing between the S 3p and Zn 3d basis functions, 
which indicates weak covalency in the Zn 3d (A^— 1) states. Photoemission studies by Langer 
and Vesely22, also indicate relatively small covalency of the Zn-S bond. A ROHF treatment 
of the type described above provides then a good description for the Zn 3d and S 3p (A— 1) 
states. 

2. Short-range relaxation and polarization effects on the valence S 3s, 3p band states 

Relaxation and polarization effects in the nearest neighbourhood of a sulphur 3s or 3p 
hole are also assessed by means of ROHF calculations on a [Si3Zn28] cluster. The cluster 
[Si3Zn28] is designed in perfect analogy with the [Zni3S28] cluster, with an active region 
[SZn4] containing a central S ion, where the hole resides, and its four NN Zn ions. The wave 
functions \^^~^a) ^'^^ ^ 3s, 3p (N — V) states are constructed by allowing for all valence 
and core orbitals within the Ca region to relax and polarize in response to the extra charge 
attached to the S 3p or 3s WO. 

The corrections to the frozen orbital (Koopmans') ionization energies due to on-site or- 
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bital relaxation turn out to be substantial, 0.85-1.30 eV, for the S 3p and 3s hole states, see 
second and third columns in Table [B As expected, the net effect of the orbital relaxation is 
larger for the S 3s hole states due to the fact that all six electrons in the polarizable S 3p 
shell are readjusting. 

The orbital relaxation and polarization effects at the NN Zn sites give rise to an additional 
moderate correction of about -0.12 eV, see fourth row in Tabled Such a moderate effect 
reflects the strongly localized, core-like character of the Zn WO's. 

Relaxation and polarization effects within the first sulphur coordination shell of the cen- 
tral S site are also explicitly considered. This coordination shell contains 12 NNN S sites. 
The distances between each of the twelve NNN sulphurs and the central S site are 3.829 A. 
Wave-function-based calculations for a cluster containing all twelve NNN S sites are compu- 
tationally unfeasible even with the smaller basis set described below. To evaluate the overall 
relaxation effect associated with the 12 NNN S ions, we designed a [S28Zn28] cluster contain- 
ing only three of the twelve NNN S sites. The Ca region [S'^Zn4S3] of the [S28Zn28] cluster 
consists of the Ca region [S'^Zn4] of the [Si3Zn28] cluster and those three sulphur sites situated 
in unit cells with lattice vectors Ri = (1, 1,0)|, R2 = (1,0, — l)f, and R3 = (0, 1, — l)f, see 
Fig. [71 They are positioned at equal distances with respect to each other and with respect 
to the central S'^ site, i.e., 3.829 A. The associated buffer region Cb contains the remaining 
9 NNN sulphur ions, all zinc ions in the first coordination shell of the 12 NNN S sites, and 
all sulphurs in the second coordination shell of the selected three active NNN S sites. The 
86-311G basis set described in Sec. Ill is employed for S where the polarization d basis 
function is, however, omitted. The basis set utilized for Zn is derived from the 86-4111d41G 
GTO basis set of Jaffe et alM by optimizing the exponents of the outermost Asp, 5sp and 
6sp shells. This optimization does not alter their values as compared to the 86-4111d311G 
GTO basis described in Sect. III. The exponents and contraction coefficients of the 3d and 
4d shells in the basis of Jaffe et al. appear to be optimal and undergo no variations in the 
basis set optimization. 

Orbital relaxation and polarization at only one of the three NNN S ions bring extra 
corrections of about -0.11 eV for both the S 3s and 3p hole states. This correction may be 
viewed as the relaxation-energy one-body increment associated with a single NNN S site in 
a many-body expansion for the overall relaxation contribution of the twelve NNN S ions^i. 

Basis set effects were found to be negligibly small, less than 0.005 eV, from additional 
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FIG. 7: Illustration of the [S28Zn28] cluster employed for the study of the relaxation and polarization 
effects within the first sulphur coordination shell around a valence-band S 3p, 3s hole. The active 
region of the cluster is the [S'^Zn4S3] fragment, where with the superscript is marked the central S*^ 
ion. The other three active S ions in unit cells with Ri, R2, and R3, and the ion are depicted 
as large white spheres. The four active Zn ions are drawn as large black spheres. The S ions in 
the buffer region are shown as smaller white spheres, whereas the buffer Zn atoms are drawn as 
smaller black spheres. 

calculations for the on-site and NN orbital relaxation and polarization on the [Si3Zn28] cluster 
and using either the 86-4111d41G or 86-4111d311G GTO basis sets for Zn and the 86-311G* 
basis set for S. The impact of the polarization d function for sulphur on the relaxation and 
polarization effects at the NNN S sites is also explicitly investigated by computing the one- 
body increment to the overall relaxation effect of the twelve NNN S ions with a [S2oZn28] 
cluster, the 86-4111d41G basis set for zinc and either the 86-311G* or 86-311G basis sets 
for sulphur. The cluster [S2oZn28] contains only one of the three NNN S sites included in 
the [S28Zn28] cluster. The active region [S'^Zn4Si] consists of the Ca region [S'^Zn4] of the 
[Si3Zn28] cluster and the sulphur ion situated in unit cell with lattice vector Ri = (1, 1, 0)|, 
see Fig. [71 The buffer region contains the other eleven NNN S as well as all Zn ions in the 
first coordination shell of the twelve sulphurs. The sulphur ions in the second coordination 
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shell of the NNN active site Si are also included in the Cb region. By utilizing either the 
86-31 IG* or 86-31 IG basis sets for sulphur, the orbital relaxation and polarization at the 
NNN active site Si give rise to a one-body increment of -0.12 or -0.11 eV, respectively for 
both the S 3s and 3p hole states. Hence, the basis set effect is found to be about 0.01 eV. 
Such basis set effects are incorporated in all numbers discussed below. 

Relaxation and polarization effects at a second sulphur site in the [S28Zn28] cluster give 
rise to an extra correction of -0.097 eV for the S 3p (N—T) states and -0.100 eV for the S 3s 
(A^— 1) states. Finally, when the orbitals at the third NNN S ion in the Ca region [S'^Zn4S3] 
are allowed to relax and polarize, an additional small contribution arises that amounts to 
-0.076 eV for the S 3p (N—t) and -0.091 eV for the S 3s (A^ — 1) states. Summing up 
the relaxation and polarization contributions of these three NNN S sites situated at equal 
distances of 3.829 A with respect to each other, we obtain an overall effect of -0.295 eV 
for the S 3p and -0.315 eV for the S 3s (A^ — 1) states. The overall orbital relaxation and 
polarization effect at the twelve NNN S sites is then calculated by multiplying by four the 
overall contribution of the three NNN S sites discussed above. The total orbital relaxation 
contributions are listed in Table HI fifth row, second and third columns. A comparison of 
those contributions for the S 3s and 3p (N — l) states reveals their small dependence of the 
angular momentum of the orbital hosting the hole. The interaction between the extra charge 
at the central S and the induced dipole at the NNN S sites is thus mainly a monopole-dipole 
type interaction. 

It is interesting to consider the relaxation and polarization effects at two NNN S sites 
which lie with respect to each other at a distance larger than 3.829 A. In the crystal structure 
of c-ZnS, this distance is either 5.415, 6.632 or 7.658 A . For this purpose, we constructed the 
[S26Zn28] and [S27Zn28] clusters, respectively. The active regions [S'^Zn4S2] of the two clusters 
contain the [SZn4] fragment and the NNN S ion located in the unit cell with Ri = (1, 1, 0)|. 
The second NNN S ions included in the Ca regions of the [S26Zn28] and [S27Zn28] clusters are 
situated in unit cells with R4 = (—1, 1, 0)| and R5 = (—1, 0, 1)|, respectively. The distance 
between the two NNN S ions in each of the two clusters is 5.415 A and 6.632 A , respectively. 
The buffer regions of the two clusters are built in perfect analogy with that of the [S28Zn28] 
cluster. All calculations were carried out with the basis sets used for the [S28Zn28] cluster. 

Calculations for the one-body increment discussed above exploiting either the [S26Zn28], 
[S27Zn28] or [S28Zn28] clusters show that this quantity is converged with respect to the cluster 
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size, i.e., a negligible increase by less than 0.003 eV is observed from the clusters [S26Zn28] 
and [S27Zn28] to the cluster [S28Zn28]. 

Orbital relaxation and polarization effects at a second NNN S site which is positioned at 
the distance d55=5.415 A with respect to a S site selected arbitrary to be the first S site, as 
in the cluster [S26Zn28], bring a correction of -0.105 eV for the S 3p and -0.107 eV for the S 
3s (A^— 1) states. When the second sulphur ion is situated at d5s=6.632 A, as in the cluster 
[S27Zn28], the extra contribution to the overall relaxation and polarization amounts to -0.106 
eV for the S 3p and -0.109 eV for the S 3s (A^— 1) states. These findings demonstrate that the 
orbital relaxation contributions associated with two distinct NNN S sites are approximately 
additive at distances between them larger than 3.829 A. The non-additive corrections to the 
overall relaxation of the joint orbital system of the two NNN S ions are small, about -0.01 
eV. We have not explicitly considered the orbital relaxation effect at a second NNN S site 
situated at d5s=7.658 A with respect to the first S site. The orbital relaxation contributions 
associated with two sulphurs at this distance are expected to be nearly additive, ~ -0.24 
eV. 

Overall contributions due to charge relaxation at three NNN S sites with distances be- 
tween any two of them larger than 3.829 A are not explicitly computed. Following a similar 
line of reasoning as for the charge relaxation at two NNN S sites, those contributions are 
expected to be about -0.3 eV. 

Finally, there is an additional contribution due to the relaxation of the orbital hosting 
the hole of -0.1 eV, see Table HI Summing up the on-site and short-range relaxation and 
polarization contributions, we obtain an overall correlation-induced shift of the S 3p valence 
bands of -2.30 eV. 

3. Short-range relaxation and polarization effects on the low-lying conduction Zn 4s, Ap band 
states 

In the following, we focus on the relaxation and short-range polarization effects on the 
diagonal matrix elements for the low-energy Zn 4s, 4p (A^+1) states. These electron-addition 
states are associated with configurations with an extra electron into the Zn 4s-like or Zn 
4p-like conduction-band WO's. To investigate those configurations, we employed a [Zni3S28] 
cluster with an active region [Zni3S4] as in the study of the Zn 3d semi core states. The 
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TABLE II: Correlation-induced corrections (in eV) to the diagonal Hamiltonian matrix elements 
for the conduction-band Zn 4s and Zn 4p added-electron states. Negative corrections indicate a 
downward shift of the conduction bands. 



A H„„(0) 



Zn 4s-like Zn 47>-like 



On-site orbital relaxation -0.068 -0.043 

NN orbital relaxation -0.908 -0.785 

NNN orbital relaxation -0.073 -0.073 

Loss of GS correlation +0.130 +0.077 

Long-range polarization -1.537 -1.537 

Total correlation correction -2.456 -2.361 



basis sets exploited are the 86-4111d311G basis set for Zn and the 86-3 11 G* basis set for 
S. In perfect analogy to the S 3s, 3p (A^ — 1) valence-band states, ROHF calculations were 
performed for the lowest-energy Zn 4s, 4p (A^ + 1) states keeping the singly-occupied Zn 4s 
or Zn 4p orbitals frozen^-. The on-site orbital relaxation and polarization in response to 
the extra electron hosted by either the Zn 4s or Zn 4p WO's amount to about -0.05 eV, 
see Table [Tll first row. This small correction to the diagonal matrix elements for the Zn 4s, 
4p (A^+1) states is consistent with the core-like character of the valence Zn 3s, 3p, and 3d 
orbitals, and the small polarizability of the Zn^"*" ion. Relaxation and polarization effects 
at the NN sulphur ions contribute to a downward shift of about 0.8-0.9 eV for the Zn 4s 
and Ap conduction bands, see Table HIl Orbital relaxation at the NNN Zn sites brings an 
additional small downward shift of about 0.07 eV of the center of gravity of the bands. 

In contrast to the S 3s, 3p and Zn 3d (N—t) states, where the WO's hosting the hole did not 
delocalize during the SCF optimization, for the Zn 4s and 4p (A^+1) states, the localization 
of the extra electron in the Zn 4s-like or Zn Ap-like orbitals could not be maintained even by 
freezing higher-energy s and p virtual orbitals at the nearest S and Zn sites. This particular 
effect is therefore neglected for the conduction-band states. 
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4- Long-range polarization effects 



An important contribution to the shift of the center of gravity of the energy bands and 
hence to the reduction of the HF gap comes from long-range polarization effects. To obtain 
an estimate of the long-range polarization effects, we compute the classical polarization 
energy of a dielectric medium outside a sphere of radius R beyond which the dielectric 
response of the crystal reaches its asymptotic value eg^: 



AE{R) 



2en R' 



The dielectric continuum approximation is a well suited approximation for the evaluation of 
long-range polarization effects beyond a given R, since the interaction of the extra charge 
with the surroundings beyond the radius R has predominantly electrostatic character, see. 



e.g., Refs. l6llJ62l . We make use of the static dielectric constant eo=8.34 deduced by Hattori 
et air- from optical data. Corrections of few eV were found to the diagonal matrix elements 
for the valence- and conduction-band states. Those corrections are included in Tables [T] and 
ITT] as well. As discussed above, relaxation and short-range polarization effects within the 
sphere of radius R were computed by ab initio methods. 

The cuttoff radii R for the Zn 3d {N — 1l) and 4s, 4p (A^+1) excited states are each 
obtained as the average of the radii of the first and third sulphur coordination shells around 
the central Zn site. Likewise, the cuttoff radii for the S 3s, 3p {N — Tj states are obtained 
as the average of the radii of the first and third zinc coordination shells around the central 
S site. Because of the equivalent structure of the coordination shells around the Zn and S 
sites, the corresponding values for the cuttoff radii R for the valence-band hole states and 
for the conduction-band electron-addition states are the same, 4.123 A. Thus, at this level 
of approximation the long-range polarization correction AH^^^^^^^{0) for each of the two 
types of states is -1.54 eV. 

The long-range polarization contributions caused by an extra charge in the conduction or 
valence bands of an insulator may be taken beyond a certain radius R to be symmetric, since 
the polarizing charges have the same absolute value. This does not hold for the short-range 
relaxation and polarization contributions (see Tables H] and HIl) because at short distances 
the associated charge distribution and polarization can not be approximated by that of a 
dielectric continuum. 
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5. Differential correlation effects: loss of ground-state correlations 

In addition to charge relaxation and polarization, a relevant contribution to the 
correlation-induced corrections to the diagonal matrix elements also arises from differen- 
tial correlation effects. An important fraction corresponds here to the loss of ground-state 
correlations. To study such effects, we construct the correlated wave functions I^I/R^mo-) 
and l^&R^Jo-) SDCI calculations starting from the reference wave functions I^R^mo-) 
I^Rrncr)) respectively. The clusters [Si3Zn28] and [Zni3S28] are employed for the study of 
correlation effects beyond charge relaxation and polarization for the S 3s, 3p (A^— 1) and Zn 
3d (A— 1) or Zn 4s, 4p (A+1) states, respectively. The active region of the [Si3Zn28] cluster 
is the [SZn4] fragment, while that of the [Zni3S28] cluster consists of a [Zn'^S4Zni2] fragment. 
The basis sets utilized in the SDCI calculations are the 86-3 IIG* basis set for S and the 
86-411 ld311G basis for the Zn ions. 

To calculate the loss of ground-state correlations, we utilized either the total SDCI cor- 
relation energies or relevant fractions thereof for the A^-electron ground state and for the Zn 
3d, S 3s, 3p (A— 1) states or Zn 4s, Ap (A+1) states. For the A^-electron GS and Zn 3d (A^— 1), 
and Zn 4s, 4p (A^+1) states, the SDCI wave functions are constructed by correlating the 3s, 
3p, and 3d orbitals of the Zn^ ion in the Ca region [Zn'^S4Zni2] plus the 3s and 3p orbitals of 
the four NN S ions. For the S 3s and 3p (A— 1) and A^-electron states as obtained from the 
[Si3Zn28] cluster, the S 3s, 3p and Zn 3d orbitals in the active region [SZn4] are explicitly 
correlated. The results for the loss of ground-state correlations are reported in Tables [T] and 
im The corrections to the diagonal matrix elements of the effective Hamiltonian due to that 
effect shift the Zn 3d bands downward by about 1.5-1.6 eV. A similar correction of about 
1.0 eV was also found for the S 3p band states. The Zn 4s and 4p electron-addition states 
are shifted upwards by about 0.1 eV due to the loss of ground-state correlations. 

In order to understand the net results reported here for the loss of ground-state correla- 
tions, some technical explanations are required without which the paper is not self-contained. 
They are found in Appendix B and should enable the reader to repeat calculations of the 
type presented here. 
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6. Overall correlation-induced corrections to the band energies 

Summing up all relevant correlation-induced corrections to the diagonal matrix elements 
for the valence-band states, the average energy position of the Zn 3d semicore states relative 
to the top of the valence S 3p bands at the F point is found to be -8.91-9.01 eV, in very 
good agreement with experimental values of -8.97 and -9.03 eV extracted from valence-band 
x-ray photoemission^^. 

The overall effect of the correlation-induced corrections considered in the current study 
to the fundamental gap is a reduction of the HF value from 10.96 eV to 5.49 eV, which is 
still 1.69 eV larger than the experimental estimate of 3.83 eV^°. We do not expect that 
more sophisticated basis sets would significantly affect those quantities. 



C. Off-diagonal matrix elements of the local effective Hamiltonian 



In the following paragraphs, we consider the correlation-induced corrections to the band 
widths. To obtain these corrections, the off-diagonal matrix elements of the effective Hamil- 
tonian in Eq. (HI) are ex plicitly computed. We apply thereby the computational approach 
adopted in Refs. 



61 



62. 



The effective Hamiltonian matrix elements between orthogonal frozen-orbital ROHF wave 
functions I^R^mo-) ^'^^ l^om'lr) having the hole at distinct sites are explicitly computed. 
These off-diagonal matrix elements H^J^^, = {^^^J^JH\^q^,]^) are directly related to the 
hopping terms of an orthogonal tight-binding representation. The same considerations hold 
for the wave functions |$RtL) and matrix elements H^^,„„,. 

Relaxation and polarization effects in the nearby surroundings of the added hole or elec- 
tron are explicitly incorporated by means of separate SCF optimizations for each of the (N—T) 
and (A^ + 1) electron states. The correlated wave functions |$R7mcr)> I'^om/a) I'^'r^V), 
l^on'a) thereby constructed. These wave functions are each expressed in terms of its own 
orbital set and are thus mutually non-orthogonal. The calculation of the Hamiltonian and 
overlap matrix elements between such nonorthogonal wave functions is performed utilizing 
an approach that has been implemented in molpro by Mitrushchenkov and Werner and 
rests on non-unitary transformations of the initial set of nonorthogonal orbitals to biorthog- 
onal sets^. 
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To construct the correlated energy bands e^^a and e^^a, one can opt for diagonalizing 
a k-dependent matrix like that in Eq. ([1]), which contains the Hamiltonian and overlap 
matrix elements between the non-orthogonal, correlated wave functions or l^^f^). 

An alternative route is based on deriving from the initial inter-site Hamiltonian and overlap 
matrix elements a set of effective hopping integrals for an orthogonal tight-binding-like 
formulation. In the latter approach, the effective hopping terms can be directly compared 
with the HF hoppings which offers a better insight into the effect of correlations on the 
electronic band structure. 

For mutually non-orthogonal electron-removal states with different binding energies, i.e., 
Ho,m'm' 7^ H^i,mm and Sn^^rnm' 7^ 0, the hoppiug terms read^ 



- s, 



Ho ,m'm' 



1 J (-^Rj,mm -^^0,m'm' 



\2 



9 1 1 — '?2 

I '-'Hi, mm' 

An analogous expression holds for tn„/(R/). AE denotes here the energy separation between 
the two eigenstates of the 2x2 non-orthogonal CI (NOCI) secular problem^. Equivalently, 
the effective hopping term can also be written as: 

When ($1-L|$^;1)=0. 

For computing NN and NNN effective hopping terms associated with the S 3s, 3p valence 
bands, two different clusters were designed, [S2oZn7] and [S26Znio]. The active regions Ca 
of these clusters are the [S2Zn7] and [SaZuio] fragments, respectively. The [S2Zn7] fragment 
consists of two S sites, denoted 1 and 2 in Fig. [HI and their NN zinc ions. Likewise, the 
[SaZuio] active region contains the S sites 1 and 2 plus an additional S ion denoted as 3 in 
Fig. [HI All NN Zn ions of these three S sites are included in the Ca subunit as well. The 
buffer regions of the two clusters consist of all S ions within the first coordination shell of 
the active zinc sites. The [S26Znio] cluster is sketched in Fig. [8] 
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FIG. 8: A sketch of the [S26Znio] embedded cluster employed for the computation of the NNN 
effective hoppings tmm' associated with the valence-band S 3s, 3p hole states. The Ca region of the 
cluster consists of the [SaZuio] fragment. The three S sites in the Ca region are depicted as large 
black spheres and labeled 1, 2, 3. S ions in the buffer zone Cb are shown as small black spheres. 
The ten Zn ions in the active region are drawn as large white spheres. 

The basis sets used in the calculations for the [S2oZn7] and [S26Znio] clusters are 86-311G* 
for the S species and 86-4111d311G for Zn. An overview of the valence-band NN and NNN 
effective hopping terms is provided in Table UTTl In the first column of Table UTTl the Wannier 
functions where the hole resides are given. The second column summarizes results from 
frozen-orbital CI calculations in terms of the HF WO's (FOCI). In the third column, we 
allow for the relaxation and polarization of the doubly occupied orbitals of the S and Zn 
ions within the Ca regions. Additionally, the hole orbital is also allowed to relax within the 
adjusted environment (ROH-NOCI), as described in Section II. 

The orbital relaxations are carried out separately for each electron-removal state I^R^mo-) 
and l^om'lr)^ constructing thereby wave functions I^R^mo-) I'^om'l)- "^^^ results indicate 
small changes of 0.01-0.04 eV for the NN effective hopping terms when relaxation effects are 
accounted for. The general trend is a slight increase in the NN and NNN hoppings tmm'- In 
each case, the relaxation of the orbital hosting the hole contributes less than 0.01 eV to the 
overall change of tmm'- The small increase of the effective hopping terms is attributed to 
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TABLE III: Nearest-neighbor (NN), R/ = (0,-1, l)f, and next-nearest-neighbor (NNN), R/ = 
(0, — 1, 0)a, effective hoppings for the valence-band S 3s, 3p states (in eV). Frozen-orbital CI results 
(FOCI ) are listed in the second column. NOCI matrix elements in terms of relaxed orbitals are 
provided in the third column (ROH-NOCI). Matrix elements smaller than 0.004 eV are not included 
in the table. 



Active WO's 




NN tw(R/) 




FOCI 


ROH-NOCI 


3s-3s 


0.080 


0.096 


3px-3px 


0.060 


0.065 


^Py{z)-3Py{z) 


0.288 


0.315 


^Py{z)-3Pz{y) 


0.368 


0.405 


•iPx{y)-'^Py{x) 


0.079 


0.088 


^Px{z)-^Pz{x) 


0.079 


0.088 


Active WO's 




NNN tw(R/) 




FOCI 


RO-NOCI 


3py-3py 


0.018 


0.020 


3pz-3pz 


0.006 


0.008 


^Py(z)-3Pz{y) 


0.005 


0.004 



the larger intersite overlap between the relaxed orbitals of the two wave functions |$^^^^,) 
and l^om'lr')- ^ valence orbitals within the NN region of the S site hosting the 

hole are polarized toward the positive charge at this site which leads to a larger intersite 
orbital overlap. Analogous trend is observed for the valence-band N 2s, 2p hole states in the 
more covalent semiconductor zinc-blende BN^. Because of the rather localized nature of the 
valence-band hole states, the effective hopping matrix elements decay fast with the distance 
R/. We found the NNN hopping terms to be an order of magnitude smaller than the NN 
matrix elements, see Table IIIII The small variations in the effective hopping integrals due 
to short-range relaxation effects lead to a small broadening of the valence S 3s, 3p bands. 

To determine the role of correlation effects beyond relaxation and polarization on the 
hopping matrix elements, we performed exploratory FO-SDCI calculations. As discussed in 
Section II, the single-particle excitations around the WO accommodating the hole describe 
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short-range relaxation and polarization effects, whereas the two-particle excitations account 
for correlation effects beyond relaxation and polarization. As for the case of diagonal matrix 
elements, we constructed CI wave functions with single and double excitations from the 
3s and 3p WO's of the active S sites. We found that the FO-SDCl treatment leads only 
to small variations of the NN effective valence-band hoppings of about 0.01 eV. A similar 
trend was observed for the valence-band N 2s, 2p hole states in zinc-blende BN— . The effect 
of long-range polarization on the off-diagonal matrix elements is not explicitly considered. 
This effect, however, is expected to be negligible. 

The semicore Zn 3d bands exhibit very small dispersion, as inferred from both 
experiment aU^>^ and theoretical results^'ii'i^'^. Our ab inito results show that the localized 
3d (A^— 1) states are stabilized by on-site and short-range relaxation and polarization, see Sec 
IV. B. 1. Since the local relaxation contributions prevail and the delocalization effects are 
very small, we do not consider explicitly hopping terms associated with the 3d hole states. 
Test calculations for the 3dt2 states show that those terms are indeed very small, about 0.01 
eV, and the corrections due to short-range relaxation and polarization are even an order 
of magnitude smaller. Those quantities are expected to be somewhat smaller for the 3de 
states. 

We carried out a similar analysis for the lower-lying Zn 4s and Ap electron-addtion states 
and conduction-band NN hopping integrals. We employed thereby the [Zn2oS4o] cluster with 
the Ca region [Zn2S7] designed in perfect analogy with that of the [S2oZn7] cluster. The Cb 
region contains all Zn and S ions in the first two coordination shells of the active S sites. 
We used the 86-311G basis for the S ions and the 86-4111d41G basis for the Zn ions. 

Results from both FOCI and NOCI calculations are summarized in Table IIV[ In the 
NOCI calculations, the valence orbitals of the two active Zn ions and the orbitals of their 
nearest S ions in Ca are allowed to relax in response to the extra electron in the conduction- 
band Zn 4s or 4p orbital. Compared to the valence-band hoppings, the larger extent of 
the conduction-band WO's determines larger values of the NN effective hopping terms. The 
correlation-induced corrections due to orbital relaxation and short-range polarization effects, 
though, are relatively small, i.e., about 0.01 eV. Variations of the hoppings occur in both 
directions: some matrix elements are slightly increased due to short-range relaxation and 
polarization, others are reduced. Analogous findings were reported for the effective hopping 
integrals related to the lower-lying conduction-band states in MgO^^ and BN— . Extra effects 
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TABLE IV: Nearest-neighbor (NN), R/ = (0, 1,— 1)| effective hopping matrix elements for the 
conduction-band states (in eV). FOCI results are listed in the second column. The third column 
summarizes the NOCI data. Matrix elements smaller than 0.01 eV are not included in the table. 



Active WO's 








FOCI 


RO-NOCI 


4s-4s 


0.081 


0.092 




1.048 


1.059 




0.260 


0.255 


^Px{y)-^Py(x) 


0.185 


0.169 


'iPx{z)-'iPz{x) 


0.185 


0.169 



arising from the relaxation of the orbital hosting the added electron are difficult to compute 
explicitly because the orbital could not be restrained from delocalizing over other cluster 
sites in the SCF orbital relaxation. Such an effect, however, is expected to be of the same 
order of magnitude, 0.01 eV, as that obtained for the ionized S 3s and 3p states. 

Next-nearest neighbour effective hopping terms are not explicitly considered since such 
calculations would require much too large clusters, e.g., [Zn26S5o], and are thus infeasible with 
the basis sets required for an accurate study of correlation-induced corrections to electron 
bands. In analogy with the NNN hopping integrals deduced for the B 2s, 2p conduction- 
band states in BN^, because of the slower decay of t„,„/ with the distance R/, some of the 
NNN tnn' for the Zn 4s, 4p states are expected to be of the same order of magnitude as the 

NN tnn'- 

Basis set effects were assessed by extra calculations at the frozen-orbital HF level for the 
NN hopping terms, using either the 86-311G basis set for the S ions and the 86-4111d41G 
basis for the Zn ions or the 86-311G* basis set for S and the 86-4111d311G basis for Zn. The 
cluster employed in these calculations is the [Zn2oS7] cluster. With the smaller basis sets, we 
found changes of about 0.01 eV in the magnitude of the NN tnn', with the largest deviation 
of 0.04 eV occurring for the NN Zn Apx-^Px hopping term. For the NNN conduction-band 
hoppings, we expect basis set effects of similar magnitude, not larger than few tens of meV. 

To summarize, the overall corrections for the intersite interactions due to the short-range 
polarization cloud around the added electron are found to be small, of the order of 0.01 eV. 
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This outcome is similar to our finding for tlie s and p-derived conduction bands in MgO^ 
and BN^. Hence, tlie widtlis of the HF valence and lower-lying Zn 4s, 4p conduction bands 
change little by incorporating short-range relaxation and polarization effects. 

Finally, we consider the effect of the correlation-induced corrections to the widths of the 
bands on the size of the band gap. We proceed by evaluating the energy ey^^ of the S 3p and 
Zn 4s, 4p band states at the F^g point. The dispersion of the px band, for example, reads 
eo,,,x=EiT^ +81^11^^^ +4tf.l^^ +0{NNN) + ...-^^ where is a diagonal matrix element 

and ti^i"'' = t^,l^^ = ti^z^'' by symmetry. Corrections of about 0.03 eV for t^x,x^^ and 0.005 
eV for t^x,x^^ (see Table [In|) lead to an upward shift of the S 3p bands by 0.196 eV at the F 
point. The band gap is reduced then by the same amount. 

We have not computed explicitly similar corrections for the Zn 3d states. Such corrections 
are expected to shift the Zn 3d bands upwards by little, such that the average energy position 
Ed of the Zn 3d states relative to the top of the valence S 3p bands at the F point will remain 
about -9.0 eV. 

The correlation-induced corrections to the Zn 4p-4p and Zn 4s-4s intersite matrix ele- 
ments are relevant at the F^ symmetry point. The energy of the 4s band at F^ is given by 
eki/,s=-E'i°i"''' + 12ti^i°^-|-...^^. For symmetry reasons, ti|l°^=ti';'i"'^\ Corrections of 0.011 eV for 
ti^s^^'' lead to a shift to lower energy of the Zn 4s band at F^ and a further reduction of the 
gap by 0.132 eV. The Zn 4p bands are also shifted downwards by about 0.08 eV. 

The correlation-induced corrections to the off-diagonal conduction-band matrix elements 
further reduce thus the size of the gap to 5.16 eV. This value is still off by 1.36 eV as 
compared to the experiment. We have excluded from the current investigation differential 
correlation effects beyond the loss of ground-state correlations. The overall contribution 
of those effects to the correlation- induced corrections to the valence-band Zn 3d, S 3p and 
conduction-band Zn 4s, 4p states is expected to lead to a further reduction of the gap. 
Studies of correlation effects beyond the loss of ground-state correlations are left for future 
considerations. Further, we have not explicitly considered the S 4s and 4p conduction-band 
states. For those states, delocalization effects are expected to dominate over local relaxation 
and polarization. As discussed in Section IV. A, the WO associated with the dispersive 
conduction band of predominantly S 4s, Ap character which is entangled with the lower- 
lying Zn 4s, Ap conduction bands is spread over four S ions and hence, very large clusters 
are needed to describe such a (A^+1) state. 
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V. CONCLUSIONS 



We have presented a detailed study of the valence- and low-lying conduction bands of zinc- 
blende c-ZnS. We have put thereby special emphasis on the "semi core" Zn 3d band states. 
The starting point in the study was a HF bandstructure. Correlation-induced corrections to 
it were computed based on a local Hamiltonian approach and a quasiparticle approximation. 
The aim of our investigation was to demonstrate that wave-function-based methods are able 
to deal in an ab initio manner with the topic of band structures without compromising on 
the infinite range of the correlation holes. This is important here, since an added electron 
(conduction band) or hole (valence band) generates a long-range polarization cloud around 
itself. Long-range polarization effects, which add to short-range relaxation and polarization, 
contribute to the reduction of the HF gap by about 1.5 eV, and therefore are not negligible. 

The new aspect in this study as compared with previous calculations on MgO and BN is 
the inclusion in the treatment of the aforementioned 3d semicore states with which density 
functional based calculations have problems. Overall correlation-induced corrections to the 
valence-band energies place the Zn 3d semicore states at about -9.0 eV below the top of the 
S 3p valence bands, in agreement with x-ray photoemission spectroscopy^^. Summing up 
all correlation-induced corrections, the size of the band gap reduces to 5.5 eV as compared 
with the HF gap of 11.0 eV. An additional, though small, reduction of the gap by 0.3 eV 
originates from correlation-induced variations in the widths of the HF bands. Correlation 
effects beyond the computed loss of ground-state correlation, when an electron or hole is 
added, are expected to decrease further the value for the gap and bring it closer to the 
experimental estimate of 3.8 eV. 

In summary, we may state that wave-function-based calculations of energy bands as per- 
formed in the present study provide a promising approach, which is free of any uncontrolled 
approximations and arbitrariness. They deserve continuing attention. 

Appendix A: 86-4111d311G basis sets for Zn. 



34 



TABLE V: Exponents (in a.u.) and coefficients of the Zn 86-4111d311G GTO basis set. 











Coefficient 


Shell 


Type 


Exponent 


s 


p 


4 




1.68680 


1.0 


1.0 


5 


sp 


0.62679 


1.0 


1.0 


6 


sp 


0.15333 


1.0 


1.0 










Coefficient 


Shell 


Tvnp 


Exponent 




d 


3 


d 


56.079664803 


029598233794 






15.748259098 


0.15879184876 






5.3098927689 


0.37984178082 


4 


d 


1.7734652510 


1.0 


5 


d 


0.51963233927 


1.0 



Appendix B: Detailed account of the loss of ground-state correlations. 

a. Zn 3(i (N — 1) states. As discussed in Sec. IV. B. 5, the SDCI wave functions for 
the A^-electron GS and Zn 3c/t2 or "idf, hole states are constructed by correlating the 3s, 3p, 
and 3(/ orbitals of the Zn^ site in the Ca region [Zn^S4Zni2] plus the 3s and 3p orbitals of 
the four NN S ions. The reference wave functions I^R^mo-) each expressed in terms of 
individually optimized orbital sets. All singlet and triplet electron pairs formed by electrons 
from the doubly occupied Zn 3s, 3p, ?>d and S 3s, 3p orbitals and from the singly occupied 
Zn 3(i(2 or 3(ie orbital are included in the correlation treatment. 

The SDCI correlation energy for the A^-electron state is found to be -22.56 eV. This 
correlation energy is partitioned into contributions of excited internal, singly external, and 
doubly external configurations. The notation for the different configurations follows that 



by Werner and Knowles in Ref. 169|. The excited singly external configurations arise from 
one-particle and semi-internal two-particle excitations from occupied to virtual (external) 
orbitals. The semi-internal two-particle excitations are two-particle excitations under the 
constraint that exactly one electron is excited from a singly or doubly occupied orbital to an 
external orbital. The excited doubly external configurations result from external two-particle 
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excitations where both electrons are promoted to virtual orbitals. The external two-particle 
excitations contribute with -22.53 eV to the total correlation energy for the A^-electron 
state. The very small contribution of -0.03 eV arising from singly external configurations is 
associated with the self-consistency of the cluster HF wave function. The coefficient of the 
reference RHF wave function in the SDCI calculation is 0.92. 

The SDCI correlation energies for the Zn 3dt2 and 3de electron-removal states amount to - 
21.54 eV and -21.43 eV, respectively. The coefficients of the reference ROHF wave functions 
in the SDCI calculations are 0.92. In analogy to the A^-electron state, the major contributions 
result from external two-particle excitations, -20.93 and -21.07 eV, respectively. For the Zn 
3d {N—1l) states, the contribution due to singly external configurations is moderate, -0.26 
eV. With the present choice of the cluster, the correlation contribution of excited internal 
configurations is about -0.09 eV for the 3de [N — Vj states and about -0.34 eV for the 3dt2 
hole states. The excited internal configurations with the Zn 3dt2 or 3de orbital doubly 
occupied are mainly related to excitations from the doubly occupied S 3s, 3p or central Zn 
3p orbitals. The smaller magnitude of this contribution for the 3de (A^ — 1) states is due to 
the orientation of the 3de orbitals with respect to the charge distributions related to the 
bonding and non-bonding electron pairs at the four S sites in [ZnS4] . 

In the framework of the quasiparticle picture and frozen local hole approximation^-, the 
occupation of the spin-orbital hosting the hole should be kept frozen in all configurational 
state functions (CSF's) in the SDCI expansions. We work, however, with spatial orbital 
rather than spin-orbital sets, and the configuration selection in molpro is done for spatial 
orbitals. Therefore, the SDCI wave functions constructed by correlating the singly occupied 
hole orbital and the doubly occupied orbitals contain CSF's where the occupancy of the 
hole orbital is no longer maintained. The contributions arising from such CSF's to the 
SDCI correlation energies are related to effects beyond the loss of ground-state correlation 
and the quasiparticle picture. In the calculation of the loss of ground-state correlations 
such contributions should be thus disregarded. These additional differential correlation 
effects may be needed, for example, for describing satellite structures in the photoionization 
spectra. 

With the present variational orbital space for the 3d (N — t) states, the occupation of 
the orbital hosting the hole remains unaltered in the excited doubly external configurations. 
Hence, the correlation contributions of those configurations are relevant in the calculation of 



36 



the loss of ground-state correlations. Contributions arising from excited internal and singly 
external configurations where the occupation of the hole orbital is altered are on the other 
hand disregarded. The energy contributions of excited internal configurations, -0.34 eV and 
-0.09 eV, are thus excluded. The correlation contribution of -0.26 eV associated with the 
singly external configurations is further analyzed in a series of SDCI calculations in which 
only certain orbitals or groups of orbitals (e.g., pairs, triples) are explicitly correlated. We 
found that a major fraction of that contribution is obtained when only the Zn 3s, 3p and 
3d orbitals are explicitly correlated, -0.21 eV for the 3d (N — t) states. When the Zn 3s, 3p 
orbitals alone are correlated, the contribution of singly external configurations is found to be 
-0.11 eV. This correlation energy is, however, discarded because of the significant presence in 
the SDCI expansion of excited singly external CSF's with both 3dt2 and 3de orbitals doubly 
occupied. Correlating explicitly only the Zn 3d orbitals results into an energy contribution 
of -0.07 eV as regards the singly external configurations. This contribution is also excluded 
from the computation of the loss of ground-state correlations. Finally, correlating the relaxed 
S 3s and 3p orbitals alone leads to a small contribution of -0.03 eV, which is also disregarded 
because of the altered occupation of the hole orbital in the singly external configurations. 
The correlation contribution of -0.26 eV is thus mainly related to configurations with an 
altered occupation of the 3d hole orbitals. Therefore, the loss of ground-state correlations 
may be deduced by comparing only the correlation contributions due to excited doubly 
external configurations for the 3d (N — l) states and for the A^-electron ground state. The 
corrections to the diagonal matrix elements due to this effect shift the Zn 3dt2 and 3de bands 
downward by 1.60 and 1.46 eV, respectively (see Table H]). 

b. S 3s and 3p (N — l) states. SDCI wave functions for the A^-electron and S 3s and 
3p (N—V) states are constructed by correlating the S 3s, 3p and Zn 3d orbitals in the active 
region [SZn4] of the [Si3Zn28] cluster. The SDCI correlation energy for the A^-electron ground 
state is found to be -24.16 eV with the largest contribution of -24.12 eV arising from excited 
doubly external configurations. 

For the S 3p (N — t) states, the SDCI energy amounts to -23.98 eV, with a dominant 
contribution of -23.30 eV arising from excited doubly external configurations. The largest 
fraction of the latter contribution, -20.93 eV, is due to correlation effects within the Zn 3d 
shell, as inferred from a separate SDCI calculation where only the Zn 3d doubly occupied 
orbitals were explicitly correlated. Excited internal configurations give rise to a small corre- 
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lation contribution of -0.13 eV and imply charge transfer from the Zn 3d doubly occupied 
orbitals, delocalized over the four Zn sites, to the S 3p hole orbital. 

A moderate contribution of -0.54 eV arises from excited singly external configurations. 
Correlating the Zn 3d doubly occupied orbitals alone brings that contribution to -0.07 eV. 
The major fraction of those -0.54 eV originates thus from correlation effects within the 
S 3s and 3p orbital space. This finding is also supported by separate SDCI calculations 
where only the S 3s, 3p and core-like Zn 3p orbitals are explicitly correlated. With the 
present cluster and such a variational orbital space, the singly external configurations are 
found to arise from excitations involving only the S 3s and 3p orbitals, with an overall energy 
effect of -0.53 eV. An inspection of the corresponding SDCI expansion indicates considerable 
presence of CSF's with altered occupation of the S 3p hole orbital as compared to that in 
the reference ROHF wave function. The correlation contribution associated with the singly 
external configurations, -0.54 eV, has been thus discarded from the computation of loss of 
ground-state correlations. To deduce the loss of ground-state correlations we utilized for the 

and S 3p (A — 1) states the dominant correlation contributions due to doubly external 
configurations. We found a correction of 0.82 eV to the S 3p band states, see Table [H 

The loss of ground-state correlations for the S 3s hole state is difficult to compute because 
maintaining the hole in the S 3s orbital in all configurations of the SDCI wave function is 
technically impossible. We calculated the associated correction by using only the correlation 
contribution of the excited doubly external configurations in the wave function I^I/r"^^) of 
the S 3s hole state. The correlation contribution resulting from excited singly external 
configurations is largely due to semi-internal two-particle excitations which involve the S 
3s and S 3p orbitals and lead to a doubly occupied S 3s orbital. It amounts to about -2.5 
eV. With the present choice of the finite cluster and variational orbital space, the internal 
configurations contribute with -0.42 eV to the SDCI energy. They are related to Zn 3d to 
S 3s charge transfer, just as in the case of the S 3p (A — 1) states. The overall correlation 
correction obtained with the selected variational space shifts downward the center of gravity 
of the S 3s bands by 1.76 eV. 

c. Zn As, Ap (N + 1) states. Analogous considerations were applied to the low-lying 
conduction-band Zn 4s, Ap states. A [Zni3S28] cluster with the active region [Zn'^S4Zni2] 
was used for the study of the (A-l-1) states and the S 86-311G* and Zn 86-4111d311G basis 
sets. We constructed the SDCI wave functions for the A^ and (A + 1) states by correlating 
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explicitly the doubly occupied 3s, 3p, and 3d orbitals of Zif and the doubly occupied 3s and 
3p orbitals of the four NN S ions. For the (A^+1) states we started from the wave functions 
l^ntncr) where the orbitals at the central Zif ion and those at the four S ions around Zn^ were 
previously relaxed. The singly occupied Zn 4s-like or 4p-like orbital is placed into the active 
orbital space but no excitations out of this orbital are included in the correlation treatment. 
The SDCI wave function constructed in this manner incorporates correlation contributions 
of external one- and two-particle excitations involving electrons or electron pairs from the 
doubly occupied orbitals as well as contributions of internal excitations and semi-internal 
two-particle excitations for which an electron is promoted from a doubly occupied orbital to 
the Zn 4s or Ap orbital. 

To evaluate the loss of ground-state correlations, we used the total SDCI correlation 
energies obtained for the and (A^-|-l) states. The corrections to the diagonal matrix 
elements Aif^^*'^(0) due to the loss of ground-state correlations are about the same for the 
Zn 4s and 4p electron-addition states, 0.1 eV. 

SDCI wave functions constructed by including in the correlation treatment not only the 
doubly occupied orbitals but also the singly occupied Zn 4s-like or 4p-like orbital account for 
correlation effects beyond the loss of ground-state correlations. Such correlation effects are 
responsible, for example, for the presence of satellite structures in the inverse photoionization 
spectra. The inclusion of the singly occupied Zn 4s-like orbital in the correlation calculation 
results, for example, into an extra differential contribution of -0.97 eV to the afore- discussed 
SDCI energy of the Zn 4s (A^-|-l) state. This extra contribution is mainly due to excited 
doubly external configurations, about -0.70 eV, but singly external configurations involving 
the single electron give rise to a relevant fraction of about -0.27 eV. We note that the 
virtual orbitals to which the extra electron is excited in the singly external configurations 
with largest CI coefficients have much weight at the four NN S ions. This extra differential 
contribution of -0.97 eV stabilizes the (A^-|-l) electron state with respect to the A^-electron 
state by -0.84 eV. It includes, however, a fraction due to a delocalization SCF effect since 
some of the projected virtual orbitals in the variational orbital space to which the Zn 4s 
electron is excited have also weight at the twelve Zn sites in the [Zn'^S4Zni2] region. We 
have not computed explicitly the magnitude of the SCF delocalization effect. Our results 
for the NN Zn 4s-4s effective hopping matrix elements at the frozen-orbital HF level (see 
Table IIVI) indicate that the delocalization of the extra 4s electron over the twelve NN Zn 
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sites is associated with a stabilization energy of about 12tS^^=-0.972 eV. 
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